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Abstract 

In this Letter we develop analytical formulations to describe sound scattering in lossless medium 
due to 2D circular wave incident on an acoustic cloak. A perfect acoustic cloak is reflectionless 
and can completely hide the cloaked object from any sound waves. However, the realization of a 
perfect acoustic cloak is difficult. Compared to plane wave, our analytic calculations show that 
circular wave from an annular line source generates distinct scattering patterns from an imperfect 
cloak design. Large modification in reflection directivities can be observed if the focal point of 
the incident wavefront is slightly customized. Hence, our work might find applications in acoustic 
cloak detection, which should have significant impact on cloak design and defense. 

PACS number: 41.20.Jb, 42.25.Fx 
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A cloak bends wave fields in desired directions to shield any interior object from detection 
l|-|4( . The theory behind cloak design is a conformal map that transforms a physical region 
with an interior cloaked hole to a virtual domain that is mathematically simply connected. 
This interesting principle is firstly proposed by Pendry et al. [5] and Schurig et al. 6| 
for electromagnetic wave cloaking. Cummer and Schurig presented a 2D acoustic cloaking 
design based on the isomorphism between acoustic equations and Maxwells equations 
Chen and Chan elegantly extended the transformation based acoustic cloak design to 3D 
spherical cloak cases Cummer et al. reported the same 3D design through a different 
route that relies on a spherical harmonic scattering analysis [7|. Acoustic cloak designs are 
generally confirmed via numerical simulations with a plane, progressive, harmonic wave. 
Some experimental demonstrations of acoustic cloak have been conducted for linear surface 
waves (at 10 Hz) 8], audible sound waves (at 1 to 3 kHz) and ultrasound waves (at 52 
to 60 kHz) 0. 

In this Letter, we develop an analytical solution that describes the interaction between 
an acoustic cloak and circular waves radiating from an annular line source. Our results 



demonstrate the existence of peculiar scattering patterns if metamaterial properties [11] 
of an acoustic cloak slightly differ from idealized ones. In addition, distinctively different 
scattering patterns can be found from circular cloak shell simply by slightly adjusting the 
focal point of the wavefront. In contrast to normally incident plane wave, this annular line 
source setup might find applications in cloak detection. 

Our investigation is conducted in an analytical way that can clearly provide physical 
insights. Before starting the analytical derivations, we have need to recall the theoretical 
foundation behind acoustic cloak. The propagation of linear sound wave perturbations is 
governed by linear Euler equations. With an e~ lbJt time dependence, these equations have 
the following form for a cloaked region, 

iujp = kkqV ■ v, (1) 

Vp = iuppov, (2) 

where p is sound pressure; v is the associated particle velocity; k is anisotropic fluid bulk 
modulus relative to k ; p is anisotropic cloak density relative to fluid density p ; p and k 
are normalized to unity. For simplicity, the discussion in the following is focused on circular 
cloaked shell in 2D polar (r — 9) coordinates (as shown in Fig. 1). The outer shell radius 



is a and inner radius is b, which is nil in the transformed virtual region for a perfect 2D 
acoustic cloak. In other words, the transformation for an ideal cloak maps the physical 
region (b < r < a) to the virtual region (0 < r < a). In the physical region, the cloak 
density is a second order tenser, that is, p = diag(p r , pg). The relative anisotropic density 
and bulk modulus are ^] 

i / \ 2 

b a \ r 



Pr = 7>Pe = ,« = r • (3) 

r — b r \a — b J r 

Outside a cloak, p r , pg and k are unity. It is easy to see that the above design is impractical 
as p r will go infinity as r approaches b. To avoid this potential singularity, the physical region 
(b < r < a) can be map ped to a virtual region (tq < r < a, ro > 0) by a linear transformation 
/. Accordingly, Norris [12| recently developed generic material specifications, 



where /' = d//dr, and d = 2 for 2D cases; d = 3 for 3D cases. 

From Eqs. (Jl])-(j2j) it can be seen that harmonic acoustic pressure (p) inside acoustic cloak 
is governed by the following wave equation 

V ■ (p-'Vp) + -p = 0, (5) 

K 

where k is the normalized wavenumber. In 2D polar coordinates, this wave equation has the 
following form, 

9/1 dp\ 1 dp 1 d 2 p k 2 , s 

o \ p r or J rp r or r z pg oa z k 

Adopting the method of separation of variables, we let p(r,9) = R(r)Q(9) and set 0(0) as 
harmonic function, 0(0) = e ime , where m is an integer. Hence, R(r) satisfies the following 
ordinary differential equation, 

d ( r dR(r)\ fk 2 p e r 2 ,\ n , , n 
From Eq. (jlj) we have p r = rf'/f and pe = f/(rf) for 2D cases, Eq. ([7]) becomes 



d 2 R(r) 1 8R{r) 



m 



It is thus easy to find that the solution of Eq. (JS} has the form: R(r) = B m J m (kf(r)) + 
C m H m (kf(r)), where J m is the mth order Bessel function of the first kind, H m is the mth 
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order Bessel function of the third kind, and B m and C m are parameters to be determined. 
As a consequence, the sound pressure in an acoustic cloak is 

oo 

P dk (r)= [B m Jm(kf(r)) + C m H m (kf(r))}e ime . (9) 



m=— oo 



Taking into account Sommerfeld radiation condition, the series forms of sound solutions 
in other regions are 



,imd 



P inc = K rnJm{kry 

m=— oo 

oo 

pSC at = J2 A m H m (kr)e md , (10) 

m=— oo 

oo 

pi nt = D m J m (kr)e ime , 



where p mc is incident sound pressure, p scat is sound pressure scattered from an acoustic 
cloak, and p mt is sound pressure in the interior of the cloaked hollow region. 
The normal velocities can be derived using Eq. fl2]), which yields 



, oo 

v? c =- K m J' m (kr)e ime , 

j oo 

vr t = _ A m (H' m (kr))e md , 
Po 

m=— oo 



V M = TO J- (B m Ji(kf(r)) + C m H' m (kf(r)))e- e , 



(11) 



m=— oo 

oo 



vT l = - E D m J' m (kry md , 
Po ^ 

' m=—oo 

where ' stands for d/dr and f'{r)/p r = (f/r) for 2D cases. Sound pressure and normal 
velocity should be continuous at cloak interfaces, that is, r = a and r = b in the physical 
region. We can therefore have the following relations, 

K m J m (kb) + A m B m {k)S) = B m J m (kb) + C m H m (kb), 
K m J' m (kb) + A m H' m (kb) = B m J' m (kb) + C m H' m (kb), 

(12) 

B m Jm{k r o) + C m H m (krQ) — D m J m {ka)^ 
^-(B m J' m (kr ) + C m H' m (kr )) = D m J' m (ka), 
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which yield 
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rn 



K, 



rn 



(13) 



D 



K, 



m 



E 2 E 6 — E 3 E 5 



where the following formulas are used, 



Ei = J m ( y kr ),E 2 = H m (kr ),E 3 = J m (ka), 

Ei = ^J^kr ),E 2 = ^H' m (kr ),E 3 = J' m (ka), 
a a 

_ E 2 E 6 — E 3 E 5 
E\E 5 — E 2 E^ 



(14) 



A 2D plane, progressive, harmonic wave can be described by p mc (x) = e tkx , which equals 
e ikrcose ^ faofc j S) Y^ = -oo i m J m (kr)e irn9 . Hence, the corresponding K m in Eq. (JTUj) is i m for 
plane wave case. Figure 1 shows a plane wave sound pressure field calculated from above 
series solutions. The 2D plane wave is incident from the left onto an acoustic cloak with 
a = 1 and b = 0.5. The normalized wavenumber is 10. For a perfect cloaking design, 
where / : (b < r < a) — > (0 < r < a), no reflection can be observed in Fig. 1(a). Since 
a perfect cloak design is difficult to implement, / : b 1— > r , r > is adopted in practical 
implementations. For example, r = 0.02 is represented in Fig. 1(b) by the dashed circle. 
However, scattering due to imperfect cloak design is now clearly visible in Fig. 1(b). 

It is of interest to study the scattering patterns if a circular sound wave is incident on 
an imperfect acoustic cloak. Figure 2 shows the setup of the problem. The circular wave 
incident on a cloak is generated by an annular line sound source, which presumably consists 
of numerous point sources obeying a uniform distribution. The radius of the annular line 
source is R, which is set to 3 in the following demonstration. The distance between the 
origin of the cloaked hollow region and the expected focal point of wavefront is L. The 2D 
annular line source gives 



that is, K m = J m (kL) for the circular wave case. Without the presence of an cloak, the origin 
of the line source will be the focal point of wavefronts. In the following demonstrations, we 
slightly move the focal point along the focal path in Fig. Practical implementation can be 
realized by marrying the concept of beamforming (in array signal processing) [13I . [yj] with 
harmonic analysis introduced in this Letter. 
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In our demonstration, Eqs. (19|)-( TI4|) can be used to calculate sound propagation and 
reflection. Figure |3] shows sound reflections due to circular waves with L = 0.1 and L = 0.3, 
respectively. It can be seen that sound reflections are symmetric with respect to the x-axis. 
The resultant scattering patterns are quite different than plane wave cases. In addition, 
distinctive scattering patterns appear if the focal point of circular wavefront is modified. 
When L = 0.1, the directivity of sound reflections are approximately in the x direction. 
When L = 0.3, sound reflections move around the y direction. For imperfect cloak cases 
in Fig. 1(b), it can be seen that almost no sound energy reaches the interior region. In 
contrast, it appears that a small fraction of sound energy arrives the cloaked hollow region. 
More details of sound transmission and reflections can be found in the online supplementary 
materials (including animations of plane wave and circular wave cases). It is for this reason 
that our work might find applications in imperfect acoustic cloak detection. 

In summary, the method proposed in this Letter permits a purely analytical design of 
acoustic cloak detection strategy. Using 2D harmonic analysis, our study demonstrates 
unique scattering pattern due to incident circular waves generated by an annular line source. 
We have also developed 3D formulations, which are omitted here for brevity. In this work 
we conducted 2D simulations using the proposed analytical formulations at various r and 
L and similar conclusions can be drawn. It is worthwhile to note that the series solution 
works to any wavelength condition. Hence, we conclude that an annular line source might be 
useful in detecting 15[ an enclosed imperfect acoustic cloak, which should have a significant 
impact on acoustic cloak design and defense. 

The preparation of this Letter is partially supported by the NSF Grant of China (Grant 
Nos. 11172007 and 11110072) and SRF for ROCS, SEM. 



* Electronic address: huangxun@pki £edu.cn| 
[1] S. A. Cummer, B.-I. Popa, D. Schurig, D. R. Smith, J. Pendry, M. Rahm, and A. Starr, Phys. 

Rev. Lett., 100, 024301 (2008). 
[2] H. Chen and C. T. Chan, Appl. Phys. Lett., 91, 183518 (2007). 
[3] S. A. Cummer and D. Schurig, New J. Phys., 9, 1 (2007). 
[4] U. Hetmaniuk and H. Liu, SIAM J. Appl. Math., 70, 2996 (2008). 



6 



[5] J. B. Pendry, D. Schurig, and D. R. Smith, Science, 312, 1780 (2006). 

[6] D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. 

Smith, Science, 314, 977 (2006). 
[7] H. Chen, B.-I. Wu, B. Zhang, and J. A. Kong, Phys. Rev. Lett., 99, 063903 (2007). 
[8] M. Farhat, S. Enoch, S. Guenneau, and A. B. Movchan, Phys. Rev. Lett., 101, 134501 (2008). 
[9] B.-I. Popa, L. Zigoneanu, and S. A. Cummer, Phys. Rev. Lett., 106, 253901 (2011). 
[10] S. Zhang, C. Xia, and N. Fang, Phys. Rev. Lett., 106, 024301 (2011). 

[11] C. M. Park, J. J. Park, S. H. Lee, Y. M. Seo, C. K. Kim, and S. H. Lee, Phys. Rev. Lett., 

107, 194301 (2011). 
[12] A. N. Norris, Proc. R. Soc. A, 464 (2008). 
[13] L. Bai and X. Huang, J. Acoust. Soc. Am., 130, 3803 (2011). 

[14] X. Huang, L. Bai, I. Vinogradov, and E. Peers, J. Acoust. Soc. Am., 131, 2152 (2012). 
[15] B. Zhang, T. Chan, and B.-I. Wu, Phys. Rev. Lett., 104, 233903 (2010). 



7 




FIG. 1: The real part of sound pressure field in the physical region. The plane wave is incident 
from the left, (a) Perfect cloak case, material properties by Eq. ([3]). (b) Imperfect cloak, material 
properties by Eq. (JIJ), / maps b = 0.5 to ro = 0.02. 
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FIG. 3: Sound scattering due to circular waves incident on an imperfect acoustic cloak by Eq. @, 
/ maps b = 0.5 to ro = 0.02. (a) L = 0.1. (b) L = 0.3, where the black solid dot stands for the 
expected focal point. 
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